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ABSTRACT 
Spectral  theory  for  linear  ideal  MHD  is  laid  on  a  firm  foundation 
by  defining  appropriate  function  spaces  for  the  operators  associated 
with  both  the  first  and  second  order  (in  time  and  space)  partial 
differential  operators.  Thus,  it  is  rigorously  established  that  a 
self-adjoint  extension  of  F(5)  exists.  It  is  shown  that  the  operator  L 
associated  with  the  first  order  formulation  satisfies  the  conditions  of 
the  Hille-Yosida  theorem.  A  foundation  is  laid  thereby  within  which 
the  domains  associated  with  the  first  and  second  order  formulations  can 
be  compared.  This  allows  future  work  in  a  rigorous  setting  that  will 
clarify  the  differences  (in  the  two  formulations)  between  the  structure 
of  the  generalized  eigenspaces  corresponding  to  the  marginal  point  of 
the  spectrum  o)  =  0. 
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1.  INTRODUCTION 

The  purpose  of  this  report  is  twofold:  1)  To  offer  to  the  applied 
mathematician  and  theoretically  inclined  mathematical  physicist  a 
self-contained  presentation  of  some  rigorous  theorems  about  the 
operators  of  linear  ideal  MHD  (linearized  Lundquist  equations)  in  a 
bounded  toroidal  domain,  and  2)  To  hopefully  attract  experts  on 
spectral  theory  and  hyperbolic  equations  to  pitch  in  a  hand  towards 
resolving  some  of  the  difficult  open  problems  remaining  concerning  the 
spectrum  of  the  linearized  operators  and  their  connection  with  the  time 
evolution  of  the  first  order  Lundquist  equations,  by  recasting  these 
problems  precisely  in  the  context  of  the  appropriate  functional  spaces. 

Concretely  speaking,  this  report  was  motivated  by  its  having  been 
pointed  out  [1]  that  there  exists  an  apparent  gap  in  the  literature  on 
the  linearized  ¥(0  introduced  by  Bernstein,  et  al.,  which  has  long 
been  known  to  be  symmetric,  as  to  whether  a  self-adjoint  extension  of 
the  latter  exists.  As  the  author  realized  this  result  was  an  easy 
extension  of  his  previous  work  [2],  he  set  out  to  exposit  this  in  the 
present  report.  However,  at  that  juncture  his  attention  was  drawn  to  a 
long  unresolved  problem  introduced  by  Harold  Grad,  as  to  the  relation 
of  the  infinitesimal  generator  and  spectrum  of  the  original  first  order 
Lundquist  equations  with  that  of  the  second  order  system  obtained 
therefrom  by  Bernstein,  et  al.  Thus  he  was  motivated  to  also  present 
herein  a  self-contained  proof  of  the  existence  of  a  serai-group  solution 


*This  work  was  supported  by  U.S.   Department   of   Energy   Contract 
No.   DE-AC02-76ER03077. 
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to  the  linearized  Lundquist  equations  and  to  thereby  provide  a  concrete 
and  rigorous  setting  within  which  the  serai-group  generated  by  them 
could  be  compared  to  the  evolution  equation  associated  with  F(5).  Some 
preliminary  remarks  as  to  the  comparison  of  the  domains  of  the  first 
and  second  order  operators  are  made  at  the  end  of  this  paper. 

We  now  pass  to  a  more  detailed  account  of  what  will  be  done  in 
this  report  and  stop  along  the  way  to  survey  the  known  results. 

The  Lundquist  equations,  or  ideal  MHD  equations,  constitute  a 
first  order  "quasi"  symmetric  hyperbolic  system  [3],  supplemented  by 
the  condition  V«B  =  0.  We  say  "quasi"  because  (see  Section  4c),  in 
general,  it  is  only  after  the  condition  V«B  =  0  is  imposed  that  the 
coefficient  matrices  are  symmetric.  In  the  context  of  nuclear  fusion 
research  they  have  traditionally  been  used  to  model  the  evolution  of  a 
plasma  in  a  toroidal  shell.  They  are  supplemented  and  coupled  with  a 
set  of  pre-Maxwell  equations  when  the  plasma  is  surrounded  by  a  vacuum 
region.  If  however  the  plasma  is  allowed  to  extend  to  the  boundary  one 
obtains  a  so-called  characteristic  mixed  initial  boundary  value 
problem.  For  the  nonlinear  case,  even  locally  in  time,  no  proof  of  the 
existence  of  solutions  exists  at  the  present  time,  and  has  appeared 
only  very  recently  for  the  characteristic  mixed  initial  boundary  value 
for  the  Euler  equations  of  compressible  fluid  dynamics  in  [4]. 

A  linearization  of  the  nonlinear  equations  around  an  equilibrium 
is  beset  with  many  difficulties  since  solutions  to  the  equilibrium 
equations  almost  certainly  do  not  exist  in  the  absence  of  symmetry  as 
was  convincingly  demonstrated  by  H.  Grad  in  [5].  Therefore  the 
following  account  which  deals  with  these  linearized  equations  will  only 
be  well  founded  for  asymmetric  perturbations  of  a  symmetric  equilibrium 
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and   needs,   as  does  all  the  literature  on  the  subject,  to  be  carefully 
interpreted  in  the  case  of  asymmetric  equilibria.   For  an   interesting 
discussion  of  these  points  see  [1]. 

We  would  like  to  point  out  that  for  axisymmetric  perturbations  of 
axisymmetric  equilibria  an  investigation  along  rigorous  mathematical 
lines  of  the  spectrum  of  F(5)  has  been  carried  out  in  [6].  Also  a 
comprehensive  survey  paper  of  results  on  the  spectrum  due  to  E.  Hameiri 
has  been  carried  recently  in  a  paper  by  the  latter  [8]. 

In  1956,  Bernstein,  et  al.,  derived,  by  introducing  a  Lagrangian 
displacement  ^  in  the  original  first  order  Lundquist  equations,  a 
second  order  symmetric  operator  F( 5)  whose  non-negativeness  allegedly 
provided  necessary  and  sufficient  conditions  for  stability.  It  was 
soon  noted,  however,  due  to  the  presence  of  essential  spectrum 
extending  to  and  including  the  origin  that  F( C)  ^0  would  be  associated 
with  non-exponential  growth  of  the  solution  [9].  Goedbloed  and 
Sakanaka  then  formulated  a  modified  energy  principle,  also  referred  to 
as  a-stability  principle  which  by  avoiding  difficulties  at  the  origin 
truly  does  provide  necessary  and  sufficient  conditions  for  a  solution 
of  pi  =  F('0  to  be  o-stable,  which  means  that  ll/p^lU  <  c  exp  at.  This 
a-stability  principle  was  put  on  a  firm  mathematical  foundation  in  [2]. 
There,  existence  of  solutions  "of  pE,  =  F(  5)  in  the  form  of  a 
evolutionary  variational  problem  following  the  work  of  Lions-Magenes, 
was  also  established.  A  second  order  semi-group  approach  to 
establishing  existence  for  the  same  equations  was  outlined  in  [10],  but 
all  details  have  not  appeared  as  of  this  writing. 

It  is  well  known  that  a  complete  resolution  of  an  operator  through 
its   spectrum  exists   only   for   self-adjoint   operators   through   the 
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so-called  spectral  theorem.  For  operators  that  are  merely  symmetric 
but  possessing  no  self-adjoint  extension,  the  spectral  problem  is  not 
well  posed  and  its  connection  with  the  time  evolution  of  the  solution 
of  the  associated  "generalized  wave  equation"  is  hard  to  interpret. 

In  the  first  part  of  this  report  we  will  show  how  machinery 
largely  developed  in  [2]  can  be  used  to  establish  the  existence  of  a 
self-adjoint  extension  of  ¥{£,)  +  XI  for  an  appropriate  X.  As  far  as 
the  author  can  tell  without  adding  a  multiple  of  the  identity  to  F( C) 
it  is  indeed  not  possible  to  find  a  self-adjoint  extension.  Moreover, 
due  to  the  the  boundary  conditions  on  £,,  the  appropriate  multiple  must 
be  determined  so  as  not  only  to  make  FCC)  positive  but  make  its 
associated  quadratic  form  coercive  over  a  certain  norm.  Of  course 
loosely  speaking  the  spectrum  of  F(5)  can  be  thought  of  as  merely  a 
translation  of  the  spectrum  of  F(5)  +  XI.  We  will  investigate  the 
self-adjoint  extension  of  F(5)  both  in  the  context  of  the  so-called 
diffuse  pinch  and  sharp  boundary  models. 

As  pointed  out  by  Harold  Grad  in  [11]  and  more  recently  in  [6], 
the  second  order  formulation  in  terms  of  F( 5)  and  the  original 
formulation  as  a  first  order  system  are  not  equivalent.  Indeed  the 
latter  is  the  more  fundamental  and  results  about  the  former  must  be 
interpreted  in  that  light.  For  completeness  and  as  this  problem  was  an 
important  motivation  for  writing  this  report  we  provide  a  brief 
account . 

The  original  first  order  system  can  be  written  following  the 
notation  in  [1] 

v^  =  L^v  (1) 
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v^  =  L2V  (2)  , 

where  u  is  the  fluid  velocity,  v  =  (p,B),  scalar  pressure  and  magnetic 
field.   Introducing  a  Lagrangian  displacement  through: 


A  hro,t)  =  vCr^.t) 


(3) 


We  then  find  by  integrating  that: 


vCr^.t)  =  L2[aro,t)J  +  (vCr^.O)  -  L2C(ro,0)j         (4) 


Thus  if: 


v(r^,0)  -  L25(ro,0)  =  0  ,  (5) 

we  obtain  from  (1): 

which  is  the  second  order  formulation  of  Bernstein,  et  al.  To 
illustrate  (but  not  exhaust)  the  difficulties  in  passing  from  (1,2)  to 
(6),  note  that  (3)  is  artificial.  An  important  class  of  initial 
perturbations  referred  to  in  [1]  as  the  anholonomic  ones  consist 
precisely  on  those  (u  v  )  which  are  not  accessible  from  an  equilibrium 
through  a  smooth  displacement. 

Another  approach  would  be  to  differentiate  (1)  and  substitute  from 
(2)  to  obtain: 
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Vj.^  =  L^L2V  =  Fv  (7) 


However,  see  [1],  only  initial  values  of  u.Uj.  are  relevant  for  which  we 
can  determine  v  from  the  relationship: 


■^^=0    ^  |t=0 


In  this  report  we  hope  to  make  a  step  towards  the  better 
understanding  of  the  difference  between  the  first  and  second  order 
formulations  by  defining  domains  for  both  which  we  conjecture  are  as 
close  as  possible.  This  is  so,  as  we  will  see,  because  the  domain  we 
will  introduce  for  the  second  order  formulation  is  the  domain  of  the 
square  root  of  F(C)  +  AI. 

Further  remarks  could  be  made  about  the  special  difficulties  in 
comparing  the  spectrum  at  the  origin  for  the  first  and  the  second  order 
formulations.   For  this  we  refer  the  reader  to  [1]. 

In  the  second  part  of  this  report  we  discuss  the  first  order 
system.  Frequently  this  system  is  coupled  with  pre-Maxwell  equations 
in  a  vacuum  region  surrounding  the  plasma.  If  however,  the  plasma  is 
allowed  to  extend  to  the  walls  (diffuse  pinch),  one  obtains  an  initial 
boundary  value  problem  with  characteristic  boundary  of  constant 
multiplicity.  A  general  theory  does  exist  for  such  problems  in  the 
linear  case,  see  [12]  for  an  up-to-date  account  and  recent 
developments.  However,  the  boundary  conditions  B«n  =  0  on  5f2  turns  out 
not  to  be  "admissible"  in  the  general  theory.  So  a  slight  modification 
of  this  theory  is  needed.  We  discuss  this  in  a  general  way  in  Section 
4c. 


-8- 

We  will  give  here  a  self-contained  proof  of  the  existence  of  a 
linear  semi-group  solution  to  the  characteristic  initial  boundary  value 
problem  assuming  only  the  Hille  Yosida  theorem  [13].  The  assumptions 
on  the  coefficients  we  will  make  will  not  be  the  most  general  possible. 
We  note  that  the  condition  V*B  =  0  for  the  system  of  Lundquist 
equations  is  nonstandard.  For  a  smooth  enough  solution  it  will  be 
satisfied  for  all  time  if  satisfied  initially.  However,  for  a  clear 
understanding  of  the  relationship  of  the  spectral  problem  with  the  time 
evolution  of  the  system  It  is  important  to  include  it  in  the  domain  of 
the  infinitesimal  generator. 

An  advantage  of  a  serai-group  solution  is  that  it  gives  a  framework, 
where  we  can  appeal  to  the  so-called  spectral  mapping  theorems  to 
relate  the  time  evolution  of  the  solution  through  the  spectrum  of  the 
associated  semi-group  to  the  exponentiated  spectrum  of  its 
infinitesimal  generator.  This  last  remark  is  especially  relevant  for 
the  Lundquist  equations  linearized  around  a  nonstatic  equilibrium. 
There,  too,  there  is  a  second  order  formulation  of  the  stability 
problem  [14],  however  in  that  case  it  is  not  a  self-adjoint  one  and 
this  further  complicates  an  assessment  of  the  relationship  of  the  study 
of  the  spectrum  of  the  second  order  operators  with  the  time  evolution 
of  the  first  order  system. 
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2.  EQUATIONS  AND  BOUNDARY  CONDITIONS 

Linearizing  the  Lundquist  equations  [17]  of   ideal  MHD  around   a 
static  (u  =  0)  equilibrium  that  satisfies: 


'Pq  =  V  X  Bq  X  Bq  (9) 


V-B^  =  0  ,  (10) 


yields  the  system: 


9v  ,   , 

1)     Pq    =  -Vp  +  VxBqXB  +  VxBxBq  (11) 


2)     I?  =  V  ^  (v^Bq) 


(12) 


3)  ||=  -YPoV-v  -  vVPq  (13) 

4)  V«B  =  0  (14) 

for  the  variables 


Here  we  have  denoted  equilibrium  quantities  with  the  subscript  o. 

When   the  plasma  region  2     extends  to  the  wall  (diffuse  pinch)  the 
boundary  conditions  for  the  equilibrium  equations  are: 


Bq.h  =  0    on    5Qp  =  <5f2  (15) 

p  =  c       on    6Q  (16) 
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and   for   the   first   order   system: 

Bti  =0  on  6J2  (17) 

vn  =0  on  5Q   .  (18) 

When  the  plasma  region  does  not  extend  to  the  wall  (sharp  boundary 
model),  in  the  vacuum  region,  Q^^  surrounding  the  plasma  we  introduce  a 
vector  potential  for  B,  VxA  =  B,  with  the  gauge  condition  V«A  =  0,  see 
[15]  for  details,  and  then  obtain  in  2^^    the  equations: 

VxVxA  =0  (19) 

V«A  =  0  (20) 

and  the  boundary  conditions: 

n^xA  =  nQx(uxBP)  ,  (21) 

-YPqV'V  +  bP.(Vx(vxBP)  +  vVBg)  =  B^.(VxA  +  vVB^)  ,    on  Tp     (22) 

and 


nxA  =0    on    r  (23) 

v 


For   the   unique   determination  of  A  in  a  toroidal  region  we  must  also 
impose  a  flux  condition: 


f   A-n  ds  =  Y(t)  (24) 

r 
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Frora  these   equations   we  may   derive,   by   introducing  a   Lagrangian 
displacement,  E,,    satisfying: 


|i  (r^,t)  =  v(ro,t)  (25) 


the  equations: 


^h 


Po  ^  =    V(yPoV-C  +   5-VPo)   +   VxBqXVx(5xBq)    -  BxVxVxC^xBq)    =:    F( 5)        in   Q      ,      (26) 
3t2 

VxVxA  =   0  (27) 

in   Q^    , 

V-A  =   0  (28) 

-YP0V.5  +   bP.(Vx(5xBP)    +    C-VBg)    =   B^.(VxA  +    5-VB^)    ,  (29) 

n^xA  =    (-n.5)B^    ,         on    Pp  (30) 

n^xA  =   0    ,  on    r  (31) 

o  V 

and  the  flux  condition: 

/   A«n  ds  =  a(t)  (32) 

r 

V 

As  described  in  [2,15],  these  equations  may  be  solved  in  the  case  of  an 
inhoraogeneous  flux  as  in  (32),  by  first  solving  them  with  homogeneous 
flux: 


/   A^n  ds  =  0  ,  (33) 

r 

V 


and  then  adding  to  A  a  solution  A'  of  the  system: 
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VxA'  =  0  (34) 

V«A'  =  0  (35) 


n^xA'  =  0  ,     on  Tp  U  r^  ,  (36) 

/   A''n  ds  =  a(t)  (37) 

V 


In  the  case  of  the  diffuse  pinch  the  equations  are: 


3t2 


=  F(0  ,  (38) 


5'n  =  0  ,     on  Fp  (39) 


We  now  associate  as  in  [2]  with  the  operator,  F(5),  a  bilinear 
form: 

a(5,i)  =  /    {Tp(V.5)(V.i)  +  Vx(CxB).Vx(ixB)  -  VxBx^.  Vx('ixB)  -  ( V*  ?)  (?.  Vp)  }dx 

+  /   VxA-VxA  d«x  -  j        (5.n)(C«n)n.V(p  +  y  Bq  ~  7  ^^  ]ds   (40) 
^v  ~    ^P 

p2    v2 

We  assume  n«V(p  +  1/2  B^  -  B^  )  >  0.   Here  A, A  are  defined  in  terms  of 

5,5   by   solving   the   boundary  value   problem   (27,28)  subject  to  the 
boundary  conditions  (30,31),  and  the  flux  condition  (33). 

In  the  case  of  the  diffuse  pinch,  the  last  two  terms  in  the  form 
are  absent  since  F,'n  -  5*n  =  0  and  the  unique  solution  of  the  boundary 
value  problem  defined  just  above  for  A  and  A  is  zero.  We  will 
hereafter  refer  to  the  diffuse  pinch  as  Case  1,  and  the  sharp  boundary 
as  Case  2. 
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3.  SELF-ADJOINT  EXTENSIONS 

In  this  section  we  will  establish  that  the  operator  F( 5)  defined 
by  (26)  has  a  self-adjoint  extension.  The  nature  of  this  extension  is 
different  in  an  essential  way  according  as  to  whether  we  are  in  Case  1 
(diffuse  pinch)  or  Case  2  (sharp  boundary).  In  both  cases  the 
extension  we  will  use  is  the  so  called  Friedrichs  extension. 

In  Case  1,  the  following  general  theorem  (see  Theorem  X23,  p.  177 
of  [16])  will  be  needed: 

THEOREM  3.1:  Let  A  be  a  positive  symmetric  operator  and  let 
a('l'.<j))  =  <'iJ,A(J)>  for  !{),<))  ^  D(A).  Then  a  is  a  closable  quadratic  form 
and  its  closure  a  is  the  quadratic  form  of  a  unique  self-adjoint 
operator  A.  A  is  a  positive  extension  of  A,  and  the  lower  bound  of  its 
spectrum  is  the  lower  bound  of  a.  Further,  A  is  the  only  self-adjoint 
extension  of  A  whose  domain  is  contained  in  the  form  domain  of  a. 

A  few  general  remarks  are  in  order  about  how  we  apply  this 
theorem.  First  of  all  in  our  case  A  =  F( 5)  is  not  necessarily  a 
positive  operator.  To  make  it  positive  we  must  add  a  multiple  X  of  the 
identity  to  F(5).  Although  we  then  satisfy  all  the  conditions  of  the 
theorem,  the  closure  of  the  form  domain  for 

axa,0   =  <F(0  +  XUO,K>  (41) 

may  not  be  such  that  all  its  elements  satisfy  C*n  =0  on  6Q.  As  we 
will  see  in  the  sequel,  however,  there  is  a  value  of  X  sufficiently 
large  (the  one  that  makes  the  associated  quadratic  form  coercive  over  a 
certain  norm)   for   which   the  form  domain  of  si^(E,,^)   does  satisfy  the 
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boundary  condition.   Thus  according  to  the  general  principle   expounded 
in   the   theorem  for   this  X,      the  associated   closed  s.a.   operator 
F(5)  +  XI,  does  incorporate   the  boundary   condition  E,'n   =  0   in   its 
domain. 

Furthermore  this  theorem  makes  clear  how  we  can  justify  in  a 
rigorous  way  a  vast  physics  literature  which  implicitly  assumes  that  a 
lower  bound  for  the  quadratic  form  yields  a  lower  bound  for  the 
spectrum  of  F(5),  a  result  which  is  not  true  for  all  extensions  of 
F(5)-  This  theorem  also  is  connected  with  the  author's  being  able  to 
justify  in  [2]  the  so  called  modified  energy  principle.  There, 
solutions  to  an  associated  evolutionary  variational  problem  were 
constructed  precisely  in  the  form  doman  of  a,. 

Case  2  (sharp  boundary)  is  unorthodox.  Indeed  a  priori  it  is  not 
clear  that  the  natural  domain  one  would  like  to  associte  with  F(5)  is 
not  empty.   This  natural  domain  is: 

{Ce  hH^I)    ■■    F(C)e  L2(J2);ypQV-5  +  bP[V  x  (  C>^bP)  +  S-VSg] 

=  B^.(VxA  +  C-VB^)}  ,     on  Fp  (42) 

where  A  is  defined  in  terms  of  E,  by  satisfying  (27,28),  (30,31)  and 
(j3).  In  [17]  Bernstein,  et  al.,  give  a  heuristic  argument  that  this 
domain  is  not  empty.  Essentially  they  argue  that  as  we  do  not  alter  A 
by  changing  the  tangential  component  of  E,  on  T  alone,  so  we  may  adjust 
this  tangential  component  in  order  that  (29)  is  satisfied.  The 
justification  of  this  argument,  involving  as  it  does  the  solution  of  a 
partial  differential  equation  on  a  closed  toroidal  surface  does  not 
appear  to  be  straightforward. 
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The  method  we  use  to  get  around  this  difficulty  is  the  following. 
With  no  a  priori  mention  of  the  operator  F( 5)  we  introduce  the 
quadratic  form  a(d)  +  X(5,5)  through  (40).  The  second  half  of  the 
theorem  quoted  above  then  insures  us  of  the  existence  of  a  self-adjoint 
extension,  ^\{K)  +  ^I(C)>  which  by  general  principles  is  defined  on  a 
dense  subset  of  L^(^Q  )  contained  in  the  form  domain  of  a(5.C)  +  A(C,C). 
By  the  nature  of  this  extension  we  have 

aU,h   +  HK,h   =  <F;^(0  +  \K,h  (^3) 

On  the  other  hand,  Green's  Theorem  implies: 

a(C,i)  +  X(5,i)  =  <F.5  +  X5,i>  +  /  (C-n)L(C,A)ds  =  0  ,      (44) 


where  L(5,A)  is  given  by  (29).  As  mentioned  above  the  set  of  E,' s 
satisfying  (29)  are  dense  in  L2(fi  )^  hence  in  H^f^p).  Thus  the 
associated  restrictions  to  the  boundary  (5*n),  are  clearly  dense  in 
L^iii   )    by  the  trace  theorems. 

This  then  implies  L(5,A)  =  0.  Hence  by  this  indirect  route  we 
have  established  that  there  exist  a  dense  set  of  E,  in  L^(Q  )  which  are 
in  the  domain  defined  by  (42). 

In  the  following  we  show  that  the  conditions  of  the  theorem  are 
satisfied  by  the  bilinear  form  (40),  and  that  the  form  domain 
incorporates  the  relevant  boundary  conditions.  Many  details  are 
omitted  and  can  be  found  in  [2]. 


CASE  1  (diffuse  pinch):  In  the  following  all  derivatives  are  taken 
in  the  sense  of  distributions.   The  domain  of  F~(5)  we  take  to  be: 


-16- 


D"(F)  =  {5eH2(n)l5«n  =  0  on  5^}  (45) 


We  introduce  the  form  core: 


W-  =  {?eHki^"p),  ^.n  =  0  on  Fp}  ,  (46) 


with  the  norm: 


"511^  =  f   {YP(V.C)2  +  IV  X  (CxBP)i2  +  X|C|2}dx  (47) 

w   J^^  o  ~ 


and  define: 


W  =  {^.n  =  0  on  FpK.V.CV  x  (5x6?)^  h'^i^^)}  (48) 


It  is  easily  seen  that  on  D~(F),a(?,5)  =  <F(C),0.   A  standard   density 

argument,   see  [18],  allows  us  to  show  {Cn^p),  C*n  =  0  on  T   }    is  dense 

in  W~  and  W.  Also  a(5,0  is  symmetric  on  W"  and  there  exist   constants 
Cj,C2,6(cpC2,X)  s.t. 

c^llCII^   <   a(C,C)   +  <5ll5llw   <   C2II51IW    .  (^9) 


see    [15]. 

The  important  point  to  check  is  that  a(5,0  is  a  closed  form  which 
is  equivalent  to  be  completeness  of  W  in  the  W  norm  defined  in  (47). 
This  essentially  reduces  to  checking  that  the  boundary  condition 
5«n  =  0  on  r  makes  sense  for  elements  in  W.  This  is  an  easy 
consequence  of  the  inequality  [18]: 
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"^•'"'H-l/2(r)  ^  ^i"5"2+  IIV-5II2}  (50) 

We   may   now  apply  Theorem  (3.1)  to  show  that  P'5^(  5)  has  a  self-adjoint 
extension  to  ^^{E,)   with  domain  given  by 


Q(Fx)  =  {C£  W|F(5)e  l2(J2)}  ,  (51) 


and  that  on  D(F^)  =  a(5,C)  =  <F^(5),p.   Here  we  have  defined, 

Fx(5)  =  F(0  +  XI(5)  .  (52) 


CASE  2  (sharp  boundary):  We  first  introduce  some  relevant  function 
spaces,  in  the  notation  of  [2].  Recalling  the  Hodge  decomposition 
theorem,  if  V  =  (Vj^,v2,v3)e  L^(^),  then  there  exists 

V^e   Vi,V2e   V2,V36V3    ,  (53) 

such  that 

V  =  V^  +  V2  +  V3  (54) 


or 


L2(n)  =  v^   +  V2  +  V3  ,  (55) 


and  the  decomposition  is  orthogonal,  where 
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V^  =  {V<t,l(t,G  H^(n)}  (56) 

V2  =  {Vxv|v  e  L2(rot)(n)}  (57) 

L2(rot)(f2)  =  {ve  L2(f2)|Vxv  eL2(f2)}  (58) 

V^  =  {h|Vxh  =  0,  V«h  =  0  in  S^,  nxh  =  0  on  652}  (59) 

All  derivatives  above  are  taken  in  the  sense  of  distributions. 

We  define  W  as  the  closure  in  the  W  norm  of  W~  =  Hl(f2  ),  where 


ICII^  =  /    {YP(V«5)2  +  IV  X  (5xBP)i2  +  X|5|2}dx  +  /    iVxAr 


Here  A  is  defined  in  terras  of  E,   by  solving  the  boundary  value  problem: 

VxVxA  =  0  (60) 

V.A  =  0  (61) 

nxA  =  {-n-O^l    ,        on  T^  (62) 

nxA  =  0  ,     on  r  (63) 

v 

/   A«nds  =  0  (64) 

^P 

The  fundamental  properties  of  W  are  given  in  Theorem  4.3   of   [2], 
which  we  recall  here. 

THEOREM  3.2:  W  is  a  Hilbert  space  of  vector  valued  function  5  for 
which 


(1)  V  X  (CxBP),  V.5,?e  L2(f2p) 

(2)  (?-n)B^e.  H-^/2(r  ) 
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(3)   The  boundary  value  problem  defined  by  (60-64)  for  A  has  a 
solution  for  which  Ae  L2(rot)(S^y)  '^  VoC^v^ 

For  the  proof,  see  [2].  Here  we  would  like  to  give  an  elaboration 
of  condition  (3).  However,  we  note,  as  will  become  clear  from  this 
elaboration  that  (3)  is  actually  a  refinement  of  (2). 

Consider  the  boundary  value  problem  (60-64).  In  Theorem  4.1  of 
[2]  the  inequality 

IIAIU    <c{IIVxAIU   +    llnxAll    _, /,         }  (65) 


I9    <c{  IIVxAlU   +    llnxAll       ,  /o         } 


was 


s  established  for  A  e  L'^(rot)(n) /^  ^2(^2),  where  f2  is  a  possibly 
multiply  connected  domain.  An  inspection  of  the  proof  of  this  theorem 
reveals  that  if  in  addition 

VxVxA  =  0  , 


then 


IIAIU    <   c  llnxAll      ,  /T  (66) 

This  follows  from  the  fact  that  if  g  is  defined  as  in  that  theorem  by 

VxVxg  =  f  (67) 

V.g  =  0  ,     in  J2  (68) 

nxg  =  0  ,     on  r  U  r   ,  (69) 

^  V 

/   g«nds  =  0  (70) 
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then 


/  VxA  •  Vxg  ds  =  0  ,  (71) 


and  we   refer   the   reader  to  the  proof  of  Theorem  4.1  to  see  how  this 
last  equality  implies  (66). 

A  consequence  of  (66)  is  that  a   solution   of   (60-64)   exists   if 
(n«5)B^  and  therefore  nxA  is  in  H"^''^(r  ).   Here  the  condition, 


/   A'nds  =  0  , 

r 

V 


must  be  interpreted  in  the  sense  of  duality 


<A,n>   ,  ,„      ,  ,„     =0  (72) 


Note,   to  write   (72)   we   have  used  again   the   result,   AS  L'^(S^), 

V«A  =  0  =>  A«ne^  H-l/2(r  U  r  ).       However,      the      condition 

V         V 

nxAS  H~^'2(p  N  ^^p^   =  0  on  r   does  not  suffice  to  ensure 

r  V 


VxAe  L^Cn^)  .  (73) 


Thus  what  characterizes  W  is  that  ^'n  lies  in  a  subset  of  H  (Tp)  for 
which  (73)  is  satisfied.  Whether  this  last  fact  entails  further 
regularity  for  £,   in  T  and/or  C  on  Tp  Is  an  open  problem. 

We   can  now  apply   again  Theorem  3.1  to  define  the  domain  of  the 
self-adjoint  operator  Fa(C)  by 


D(F^)  =  {5e  W|a(5,C)  =  <Fx(C),5>} 


2 
The  theorem  guarantees  that  this  domain  is  dense  in  L  (Q) 
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4.   FIRST  ORDER  LUNDQUIST  EQUATIONS:   SEMI-GROUP  SOLUTIONS 
AND  SPECTRAL  MAPPING  THEOREMS 

To  begin  this  presentation  we  make  a  change  of  variables  that  puts 

the  Lundquist  equations  in  canonical  symmetric  hyperbolic  form: 


iit  =  I  Ai  -^  u  +  Bu  (74) 


u   an  n-vector,  A.  symmetric  n-n  matrices.   We  do  this  in  the  following 
subsection. 

4a .   Quasi-Symmetric  Hyperbolic  Form  of  Equations 
Let 


P  =  (YPo)"^^^P  (75) 

V  =  (Po)^/2v  (76) 


In  the  new  variables   the  Lundquist   equations   linearized   around   an 
arbitrary  equilibrium  (static  or  nonstatic)  can  be  re-written: 


V  =  -[I^jl/^V;  -.  '^"°Y"'^'  Vp,;  +  (p,)-l/2v.BxB,  +  p-l/2vxB,xB 
Po  ^ 


^o  PQ 


B  =  Po^/^V  X  (vxB)  -  1  p-3/2  [Vp^  x  (vxB) ]  (78) 
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■YP 


Pt  =  -(-^]^/2v.v  +  i.  (YPoPo>^^^^*^Po  -  (YPoPo)"^^^^*^Po 


-  Vq-v;  -  l^-^v^-Vpo  (79) 

\  =  -Po^^^  ^'^^o  -  %*^S  (80) 

with  principal  part 

v^  =  -[_^J^/2vp  +  p-1/2vxBxBq  +  ...  (81) 

Pq 

B^  =  Po^/^Vx(;xB)  (82) 

Y  D 

P^  =  -[-^jl/2v.v  +  ...  (83) 

Pq 

St  =  -Vq-VS  (84) 

To  this  set  of  equations  we  must  add  the  condition: 

V.B  =  0  (85) 

In  this  form  with  the  help  of  a  few  vector  identities  the  symmetry 
of  the  A^  is  easily  ascertained  directly  when  Vq  =  0.  We  next  discuss 
in  a  general  way  the  relationship  of  the  semi-group  solutions  we  will 
define  with  the  spectrum  of  their  infinitesimal  generators.  The 
discussion  will  apply  whether  the  underlying  equilibrium  state  is 
static  or  non-static  (with  flow).  We  hereafter  drop  the  tilda  in  the 
perturbed  quantities. 
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4b .   Remarks  on  Spectral  Mapping  Theorems 

The  assumed  smoothness  of  Pq.Pq  and  B^  allow  us  to  define  all  the 
spatial  differential  operators  arising  in  (74)  in  the  sense  of 
distributions.   Let 


V 

B 


and  let 


Lu  =  A.  u  +  Bu  , 

1  9x. 


where  the  derivatives  are  taken  in  the  sense  of   distributions   and  A- 
and  B  are  determined  by  (77-80). 

A  natural  definition  for  the  domain  of  L  would  be 


D(L) 


u  e {L2(n)}7  :  V.B  =  0  ,  Lu  &  {l2(Q)}7|, 
v°n  =  0  on  6Q,    B«n  =  0  on  5i2 


(86) 


We  hereafter  will  use  the  shorthand  uG  L^{n)  instead  of  u  e  {L^iQ)}' . 
On  this  domain,  however,  we  conjecture  that  L  is  not  a  closed  operator. 
This  is  due  to  failure  in  general  of  the  condition 


Il4i!-^A,  >0  (87) 


Thus  we  will  instead  consider 
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L  +  AI  , 

where   A  is  determined  so  as  to  ensure  that  L  +  AI  is  a  closed  operator 

if  A  >  A  .   It  can  also  easily  be  shown  that  for  the  proper   choice   of 

Aq   all   the   domains  of  the  associated  L  +  Al  are  the  same.   Thus,  for 

the   purposes   of   this   preliminary   discussion,    without    loss    of 

generality,  we  fix  a  particular  A  >  A  . 

L  +  AI   will   be   shown   to   be   the   infinitesimal  generator  of  a 

semi-group  denoted  T_(t).   Formally: 

\ 


.(L  +  AI)t  ^  T_(^)  ^  ^At^Lt  ^  e^S.(t) 
A 


(88) 


where  T(t)  is  the   solution   to   original   system   (74).    Thus   for   a 
particular  initial  data  u   the  solution  of 


IS 


u.  +  y  A-i  ^  u  +  Bu  =  0  (89) 


u(x,0)  =  u   ,  (90) 


u(x,t)  =  e-^'^  T_(t)Uo  ,  (91) 

X 


for  any   u  e  D(L  +  Al).    We  recall,  as  is  well  known,  [13],  that  the 

set  of  A's  in  the  spectrum  of  T(t)  =  e~^  T_(t)  are  contained  In  a  strip 

A 
in   the   complex  plane,  a_  <  Re  A  <  a+.   Once  we  have  shown  that  L  +  Al 

generates   the   semi-group   T_(t),   the   spectral   mapping   theorems 

A 
[19,  p.  44-48],  will  guarantee  that  the  following  set  Inclusions  hold: 
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{gOCL  +   XDtj  (^  a[T_(t))    ,  (92) 

X 


therefore. 


.-Xt^g[a(L     +   Al)tj  c^g-Xtg^x_(t)]    ,  (93) 

X 


but, 


sup        _        e^^^  "^   ^^^"^    <    IIT_(t)ll  (94) 

ae  a(L  +   XI)  X 


Now, 


sup    (|X|  |x  e  e'''^'^a[T_(t)]}    <   e"-^*^  IIT_(  t)  II  (95) 

'  X  X 


and  therefore 


sup_    e'^'^e^^  +  ^^^^  ^  e-^^IIT_(t)ll  (96) 

(L  +  XI)  X 


Equality   holds,  for  instance,  when  T_  is  self-adjoint,  which  it  is  not 

X 

in  our  case. 

The  foregoing  inclusions  are  the  only  general  statements  that  can 
be  made  in  abstract  way  as  to  the  relation  of  the  spectrum  of  the 
infinitesimal  generator  and  the  time  evolution  of  the  associated 
serai-group.  They  show  that  we  cannot  expect  to  get  from  the  spectrum 
of  the  infinitesimal  generator  necessary  and  sufficient  conditions  for 
stability.    This   is  moreover  not  due  to  the  lack  of  regularity  of  the 
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semi-group  solution  as  the  problem  exists  even  for  analytic  serai-groups 
[20]. 

The  spectral  mapping  theorems  do  however  justify  through  (96),  the 
procedure  of  deriving  necessary  conditions  for  stability  or  put  a 
different  way  sufficient  conditions  for  instability  from  the  knowledge 
of  the  range  of  the  spectrum  of  L  +  XI. 

When  linearizing  around  a  static  equilibrium  we  can  say  a  good 
deal  more  at  least  In  a  heuristic  fashion  about  the  time  evolution  of 
the  Lundqulst  equations  by  considering  the  spectrum  of  the  self-adjoint 
second  order  operator  F(5),  see  (26),  which  arises  when  we  integrate 
the  system  (2)  so  as  to  obtain  an  equation  for  ^  alone.  Heuristically , 
we  would  expect  that  any  exponential  growth  of  the  semi-group  T(t) 
would  correspond  to  a  negative  element  in  the  spectrum  of  F(5)« 
However  there  are  difficulties  In  justifying  this  argument  as  we 
pointed  out  In  the  Introduction,  which  we  will  not  resolve  here  beyond 
the  comments  made  at  the  end  of  this  work,  as  to  the  comparison  of  the 
domains  of  F(5)  and  L.  What  we  would  like  to  stress  at  this  juncture  is 
that  when  linearizing  around  a  non-static  equlllbrlm  where  there  Is  no 
self-adjoint  second  order  formulation,  a  sharper  version  of  the 
spectral  mapping  theorems  may  prove  the  most  viable  route  towards  a 
better  understanding  of  the  relation  of  the  spectrum  of  L   (or   L  +  XI) 

to  the  time  evolution  of  the  semi-group  T(t)  [or  T_(t:)j. 

X 

We   now   turn   to   the   proof   of   the   existence   of  a  semi-group 

solution.  We  will  consider  separately  the  cases:  1)  linearization 
around  a  static  equi  librin-n,  2)  linearization  arou'id  an  equlllbriam 
with  flow,  and  3)  incompressible  plasmas.  Before  we  do  this,  we  will 
make   some   general   remarks  about  the  non-applicability  of  the  general 
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theory  of  maximal  positive  boundary  data  for  symmetric  positive  systems 
to  our  case.  These  remarks  may  be  of  interest  only  to  those  already 
familiar  with  the  theory  and  may  be  skipped  without  loss  of  continuity. 
We  consider  the  static  case  separately  as  opposed  to  a  special 
case  of  the  non-static  case  in  order  to  facilitate  future  work 
regarding  the  comparison  of  the  first  and  second  order  operators  and 
associated  questions  mentioned  previously. 

4c.   On  the  Inapplicability  of  General  Theory  for 
Symmetric  Positive  Systems 

This  section  is  intended  for  those  who  are  familiar  with  the 
theory  of  maximal  positive  boundary  data  for  symmetric  positive 
systems. 

We  recall  the  general  principle.  Suppose  we  denote  by  N(x)  the 
linear  subspace  in  which  the  vector 


must  lie  in  if  it  is  to  satisfy  the   boundary   condition  at   x.   N   is 
called  maximal  positive  if 


<A  (x)w,w>  >    0      ,        V  X  e  6^  ,    w^NCx)  ,  (97) 


where 


An  =  A^n^  ,  (98) 
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and  if  N  canot  be  enlarged  while  preserving  (97).    Note   that   clearly 
ker  ^C   N. 

Now  a   simple  calculation  shows  that  the  matrix  A   for  the  system 
(77-80)  has  the  form, 


(-n-v„) 


-[]^jl/2n.    p-l/2[(n.B^)-n(B„.) 


(-Vo-n) 


p-l/2[(B^.n)-BQ(n.)] 


-(vQ«n)  +  VgCn*) 


/ 


-v^.n 


The  above  matrix  is  so  far  completely  general.  That  is,  we  have 
not  yet  imposed  any  boundary  conditions  on  the  equilibrium  quantities 
Vq,  Pq,  and  B^.  We  have  used  the  notation,  B  (n«),  to  denote  the 
operator  submatrix  that  acts  on  w  =  (w.,wn,wo)  through. 


(99) 


(B^nJ)(w)  =  BjnJwJ  . 


(100) 


The  transposed  matrix  operator  is. 


<",. 


(101) 


which  we  denote  by. 


n(B^.) 


(102) 
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The  matrix  (99)  is  manifestly  symmetric  except  for  the  expression, 

-(v^.n)  +  Vo(n.)  ,  (103) 

whose  transpose  is, 

(-Vq.h)  +  (nCvQ.)  (104) 

We  see  that  (103)  and  (104)  are  not  generally  equal.  This  is  why  we 
refer  to  (77-80)  as  a  quasi-symmetric  hyperbolic  system.  When  we 
restrict  the  action  of  (99)  to  elements  B  satisfying  B«n  =  0  the 
symmetry  is  salvaged.  B«n  =  0  is  valid  for  B^  that  are  subjected  to 
the  condition  V»B  =  0. 

Note,  however,  that  for  static  equilibria  v  =  0,  and  the  symmetry 
holds  even  without  V«B  =  0. 

We  now  turn  to  the  discussion  of  the  matrix  (99)  when  it  operates 
on  the  boundary.  There  the  conditions  n«v  =  n»B„  =  0,  p^.  =  c  must  be 
imposed.  It  is  to  the  corresponding  simplified  matrix  and  its  operator 
kernel  that  we  will  now  refer.  Again,  even  at  the  boundary,  this 
matrix  is  only  symmetric  when  acting  on  a  domain  consisting  of  u's, 
u  =  (v,p,B),  for  which  V«B  =  0.  Even  on  this  restricted  domain  the 
conditions  (97)  of  the  general  theory  do  not  hold  however.  To  see  this 
consider  an  element  u  =  j(v,p,B,S)  which  satisfies, 

n-v  =  0  (105) 


^  (Bo-B)  +  [^]'^^V   =  0  (106) 


o 
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It   is   easily   seen   that   such  u   satisfy  uS  Ker  A    However 
although  v»n  =  v«n  =  0,  (106)  does  not  imply  that  B«n  =  0,  thus 

Ker  fi^C^   N(x)  .  (107) 

Of  course  reinstating  the  condition  V«B  =  0  and  its  "characteristic" 
equivalent  n»B  =  0, 

Ker  A^c  N(x)  , 

if  we  interpret  A^  as  acting  only  on  this  subspace.  Note  that  given 
the  change  of  variables  (75-76),  (106)  is  the  linearized  version  of  the 
condition, 

p  +  b2  =  c  ,  (108) 

which  thus  appears  to  be  a  more  natural  boundary  condition  than  the 
traditionally  imposed  one  B«n  =  0. 


4d.   Existence  of  Semi-Group  Solutions 
We  first  recall  the  definition: 

DEFINITION:  A  semi-group  T( t ) ,  0  <  t  <  «>  of  bounded   operators   is 
said  to  be  a  continuous  C^  serai-group  in  a  Banach  space  x  if: 


lim  T(t)x  =  X  ,     for  every  xQ  X 
t+O 
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These  remarks  may  be  of  interest  only  to  those  already  familiar  with 
the  theory  and  may  be  skipped  without  loss  of  continuity.  We  also 
recall  the  Hille-Yosida  theorem: 

THEOREM  4.1 :  An  operator  A:  X  -»■  X  is  the  generator  of  a  serai-group 
of  bounded  operators  T(t)  satisfying 

IIT(t)  II  <  ce'^'^  , 
if  and  only  if 

(1)  A  is  a  closed  operator  with  dense  domain. 

(2)  The  resolvent  set  p(A)  of  A  contains  (a),»). 

(3)  U(AI-L)"ll  <.   M/jX-ojl". 

Furthermore,    if  M  =    1    can   be   deduced   from 


3'  lUI-AII    <  I —  (109) 

|X-a)l 


We  will  show  that  the  conditions  of  the  theorem  are  verified  for  the 
operator,  L  +  XI  defined  by  (77-80)  on  an  appropriate  domain.  Details 
will  be  given  only  for  linearization  around  a  static  equilibrium. 

CASE  A  (static  equilibrium):  X,  the  underlying  Hilbert  space,  will 
be  defined  by; 

X  =  {v  =  (u,p,B)|u  eL2(J^),  V-B  =  0}  (110) 


-32- 


REMARK:  Note  that  L  maps  (v,B,p)  to  (-,  V  x  (v/p^^  x  B^),  -).   As 


V-[V  X  (_^  X  Bq)J  =  0  ,  (111) 

Po 


Lu  e  X,  provided  Lu^  h'^i'il) .      Clearly  the  same  holds  for 


Lu  +  XI  (112) 


X  will  be  determined  by  (77-80).   We  now  define: 


D(L_)  =  {u  =  (v,p,B)|ue  L'^(n),  L_ueL''(J^),  V.B  =  0,  vn  =  0,  B.n  =  on  5^}   (113) 
X  X 

Clearly  D(L_)ex.   Furthermore,  if  6n  is  smooth  enough,   say,   C  , 

X 
and  we  define, 


D+(L;^)  =  {u  =  (v,p,B)|vecl(^)  , 

Vp  -  VxExB^e  C^(:2),  V.B  =  0,  v'n  =  0,  B.n  =  0  on  6^} 


D'''(L_)  is  dense  in  D(L_)  in  the  graph  norm  IIL_u  II .   To  show  this  we  will 

XX  X 

first  need  to  show 


THEOREM  4.2:  L_  is  a  closed  operator  for  sufficiently  large  X. 
X 


PROOF:  We   first   show   that   on   the  domain  D"*",   the   following 
inequality  holds  for  larger  X: 
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1L_  +  Xull  >  cllull  (114) 

A 


For  ueD"*",  this  is  a  consequence  of  the  symmetry  of  the  A. 's  and  the 
boundary  conditions  u«n  =  B«n  =  0  on  5fi.  One  reasons  as  follows:  Let  X 
be  real,  consider 


<(L  +   XI)u,(L  +   XI)u>   =    llLul|2  +    X^llull^  +   2   Re    X  <Lu,u>  (115) 


Now 


<Lu,u>  =  <A^  .JL-   u  +  Bu,u>  (116) 


cJ^   u,Aj^u>  +  <Bu,u>  (117) 

3x^ 


<u,  -A_  (A.u)  +  <Bu,u>  (118) 


<u,  [^   A.  lu>  -  <u,Ai  -J—   u}  +  <Bu,u>  (119) 

9x.  ^■^  1  dX. 


Here   in   the   second   step  we  have  used  the  symmetry  of  the  A^,  in  the 
third,  the  integration  by  parts,  the  boundary  terras  dropped  out  because 


/   <Ainj^u,u>ds  =  0  (120) 


Specifically  the  boundary  terms  that  arise  are 


P(vn),  {(vxB  )  X  B}  -n  (121) 


-34- 
The  second  terra  can  be  written 

(vQ.B)BQ.n  -  (BQ.B)v«n  ,  (122) 

which  vanishes  by  virtue  of  (17)  and  (18).   It  is  noteworthy   that   the 
boundary  condition, 

B«n  =  0  , 
was  not  used  to  establish  (118).   Thus,  from  (118)  we  obtain 


2  Re  <A.  _i_  u,u>  +  2  Re  <Bu,u>|  <  c'  llull^  (123) 

8x. 


Thus  for  X  >  c'  +  X, 

ll(L  +   XI)ul|2    >    IILull^  +    Xllul|2    .  (124) 

Moreover,  an  elementary  calculation  using  the  system  (77-80) 
establishes  that  for  X^  sufficiently  large  there  exists  d  and  X  so 
that,  if  X  >  Xq 

llLul|2  +   Xllull2    s   d{ll-Vp  +   VxBxB^^||2 

+    IIV   X    (vxBq||2  +    llV«vl|2  +    Xllul|2}  (125) 

We   define 


L_   =   L  +    XI    ,  ^   >    ^o 

X 
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It  is  clear  from  the  inequality  (124)  that  the  norm  llL_u  II  is  equivalent 

to  the  graph  norm  {IIL_ull^  +  llull^}^'^  on  D"*"  and  thus  on  the   closure   of 

X 
D"*"  in   this  norm.   Thus  to  show  that  (113)  is  valid  on  all  of  D(L_)  we 

X 
need  merely  to  show  that  D"^  is  dense  in  D.   For   now,   assume  we   have 

proved   this.    We  are  then  in  a  position  to  conclude  the  proof  that  L_ 

X 
is  a  closed  operator. 

Suppose  that  (u^,L_Uj^)  is  a  Cauchy  sequence   in   the   graph  norm. 

X 
Since   X   is  complete,  there  exist  u,V  such  that  u^  ->■  u,  Vj^  =  L^Uj^  ->•  V. 

We  wish  to  show  that  u^  D(L_)  and  L_u  =  V.   With  u  =  (v,p,B) 

X       X 

V  X  (vxBq),V«v,  Vp  -  VxBxBoS  L^(^)  ,  (126) 

V.B  =  0  ,  (127) 


are  consequences  of  the  continuity  of  differentiation  is  the  sense  of 
distributions.  So  we  need  merely  to  check  that  the  boundary  conditions 
B«n  =  0,  vn  =  0  are  satisfied.   As 


Luj^  >   c'    >   c'  IIBn"  (128) 


B«n  =  0  on  6J^  is  a  consequence  of  the  inequality 


IIBII      ,  /„  <   cllBII  (129) 

H-l/2(r) 


Similarly,  from  this  inequality  applied  with  u  instead   of   B   combined 
with  the  fact  that 

IIL;^(un-Un,)  II  >   c'IIV.(Vn-Vn)ll  ,  (130) 
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we   can   conclude, 


vn  =   0    ,  on    5J2   .  (131) 


thus    ue   D(L_),    L_(u)    =   V. 

A  X 

Q.E.D. 


We   next  would   like   to   verify   condition  2   and   3'   of   the 

Hille-Yosida  theorem.    Condition  2,   that   the   resolvent   set   pC^x) 

contains   (oj,")   for   some   sufficiently   large  cj  follows  easily  from  a 

general  argument  which  we  outline  here  for  completeness.    Firstly,   as 

L_   is  a  closed  operator,  clearly  the  range  of  L_  on  D(L_)  is  closed  in 

X  XX 

X.  Next  we  note  that  for   X^   sufficiently   large   and   real   [possibly 

larger   than   X^  defined  by  (124)j,  (113)]  applies  to  the  adjoint  L^  of 

L^  to  show  that,  if  X'  >  Xq,  uS  D(L;^) 


l(L*  +   X'Dull    >   clluil  (132) 


This  shows  that  the  range  of  L^  +  X'l  is  dense.   For   if   it  were   not 
there  would  exist  VG  X,  V  ^  0,  for  which 


<(L_  +  X'I)u,V>  =  0  ,  (133) 

\ 


This  would  then  be  true  in  particular  for  VGC  (JJ).  Taking  into 
account  the  definition  of  L,  +  X'l  as  a  differential  operator  acting  in 
the  sense  of  distributions,  this  then  implies  v^  ^^^x^  and, 


(Lt  +  X'I)V  =  0  ,  (13A) 
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which   by  virtue  of  (132)  implies  u  =  0.   Lastly  we  note  that  condition 

3  of  the  Hille-Yosida  theorem  is   satisfied   by   virtue   of   inequality 

(113). 

There  remains  only  to  prove  the  density.   We  recall  the  definition 

of  D(L_), 
X 


D(L_)  =  {ug  L^{^),    -Vp  +  VxBxBq,  V  x  (vxB^), 
\ 

V«ve  L2(n),  V«B=0,  v«n  =  B«n  =  0  on  5f2} 


The  technique  in  [18]  applies.  See  Remark  [4.2]  and  pgs. 
339-341.  In  outline  one  seeks  to  find  approximating  sequences  that  are 
defined  on  an  open  set  Q' (Z  f2.  The  elements  of  this  approximating 
sequence  can  then  in  turn  be  approximated  in  Q  by  first  multiplying 
them  by  a  function  which  is  one  Q  and  zero  outside  Q' ,  and  then 
mollifying  the  result.  To  define  the  original  sequence  one  uses  the 
existence  of  an  open  covering  of  the  set  Q  and  SQ  whose  covering 
property  is  invariant  under  translations.  For  a  more  sophisticated 
approach  under  weaker  hypothesis  on  the  coefficients  see  [12].  Then 
the  translational  continuity  of  L^i^l)   does  the  rest. 

We  next  would  like  to  outline  how  to  handle  the  cases  2) 
equilibrium  with  flow  and  3)  incompressible  plasma. 

CASE  B :  Equilibrium  with  Flow  (Diffuse  Pinch) 

We  first  note  that  in  the  presence  of  flow,  the  equilibrium 
equations  become, 

Po^o-^^o  =  -  VPo  +  V  X  Bq  X  Bq  (136) 


-38- 


%-VPo  +  YPoV-Vo  =  0  (137) 

Vq.VSq  =  0  (138) 

V«Bq  =  0  (139) 


p  may  be  determined  from,  (Sq.Pq)  from  the  equation  of  state 

Po  =  %pI  (^^0) 

The  boundary  conditions  that  are  imposed  are 

B^.n  =  v^n  =  0  ,     on  r  (141) 

We   explain  now  how   the   special   structure  of   the   Lundquist 

equations  makes  possible  the  definition  of  a  closed  operator  L   on  a 

domain  that  includes  the  boundary  condition  v«n  =  0,  B«n  =  0. 

First,  we  again  define  for  an  appropriate  X 


Sue  l2(S^)  ,  L^uG.L2(il)   i 
v«n  =0   on  6n  ,  V.B  =  0  ^ 


(1^2) 
D(L,)  = 


To  see  that  this  domain  is  closed  in  the  graph  norm  we  must  show  that 
the  condition  vn  =  0  is  satisfied.  All  the  other  steps  are  similar  to 
those  in  the  case  of  static  equilibria. 
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Let, 


Lu 


r-  L,u 


L2U 


L^u 


(143) 


L/.u 


where  Lj^,L2,L3,L4  correspond  to  v,.,  B^.,  pj-   and   Sj-   respectively.    By 
definition, 


Vq 

L,u  =  V  X  [-1—   X  B_l  +  V  X  [__-  X  B] 


(144) 


Now  given  B,  note  the  identity 


V-((-^  X  Bq)  X  b)  +  '^•[(-j^  X  B)  X  i] 


=  B-V  X  [_Z_  X  Bj  -  [-^  X  Bj  .  Vxi  +  i-VP  X  [_!^  X  BJ 


[ °-   X  Bl  •  7xB 

pl/2     ^ 
^o 


(145) 


Integrating   over   Q  and  using  the  divergence  theorm  we  obtain  for 
the  L.H.S.   of  (122) 
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V  


(146) 


Po  Po 


which    upon    being    expanded    and    using    the    conditions 
B«n  =  B^.n  =  v^'n  =  0  yields, 


/   (v-n)BQ.B  ds  (147) 


Choosing     now     B  =  aB^,     with     aeH^'^(r),     such    that 
Hall  i/o(r)  >  cllall  1    ,  we  find  from  (124)  and  (122)  that, 

\j      (vn)B2  a   ds  | 

^  l^^o'l  CO  "a"}  11^2""  +  ^"'  ""  '^"  (148) 

Li 

<  cjllall    IIL2UII   +    llVa  X   BqII2    HvH   +    Hall    llv  II  j  (149) 

<  c    I  Hall    lILull   +    Hall    ,  Hull}  (150) 

h1(J2)  ' 

<  c|  Hall    ,  IILull  +    Hull}  (151) 

<  c'  HaH    .  ,„         j  IILuU    +    HuH}  (152) 

Hl/2(r)'  J 


Thus 


IvnB^ii         .  i,   c' I  IILull   +    IIuHI    .  (153) 


Now  since  B^vn  =  0  in  H  ^'^(T)  implies  vn  =  0  in  H  ^'^(T)  (given 
that,  by  assumption,  B^  #  q  on  F),  we  see  that  the  condition  vn  =  0  is 
preserved  by  this  norm. 


-Al- 


CASE  C :  Incompressible  Plasma 

Here  we  add  the  condition  V«v  =  0  in  U.   In  this  case   it   is   not 
appropriate  to  make  the  change  of  variables  (75-76).   We  use  only 


,1/2 


P  =  (ypo)''^  P  (154) 


This  transformation  puts  the  system  (77-80)  into  the  form 


(p_I)u^  +  a:.  ^  u  +  B'u  =  0  (155) 

o    '^     1  9x. 


with  symmetric  A^,  or 


u^  +  A-  -^  u  +  B'u  =  0 
t    1  3xi 


4'  =  (pol)-l  Ai  .  B^'  =  (p^I)-l  B^ 


(156) 


Thus  A£'  are  also  symmetric.   It  is  now  easy  to  prove  that 


L,  =  -A^'  --i-  +  B"  +  XI  ,  (157) 

^  1   8x^ 


generates  a  serai-group.  The  boundary  conditions  vn  =  B«n  =  0  on  6J^ 
now  both  are  preserved  in  the  graph  normal  due  to  the  inequality  (50) 
and  the  fact  that  V«v  =  V«B  =  0. 
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5.  COMPARISON  OF  THE  DOMAINS 

We  here  compare  the  domains  of  the  operators  corresponding  to  the 
second  and  first  order  formulations  for  the  case  of  the  diffuse  pinch. 

As  our  notion  of  solution  to  the  second  order  equations,  as 
discussed  in  [2]  is  to  associate  with  thera  a  second  order  evolutionary 
variational  problem  following  the  work  of  Lions-Magenes.  The  solution 
is  thus  defined  in  a  space  of  functions  requiring  L^  restrictions  on 
only  the  first  order  derivatives.  This  is  what  accounts  for  the 
similarity  between  the  domains  (48)  for  the  second  order  and  (135)  for 
the  first  order  (symmetric  hyperbolic)  formulations.  As  is  clear  from 
these  definitions  the  space  defined  for  the  first  order  theory  is  such 
that  its  first  component  vlies  in  the  space  W  defined  for  the  second 
order  theory.  If,  instead  of  v,  we  require  as  discussed  in  the 
introduction,  that  5&W,  and  define  p,  B  which  were  denoted  in  (2)  by 
V  by 


p  =  -C'VPq  -  YPo^*?  +  gl(x) 

(158) 
B  =  V  X  (CxBq)  +  g2(x) 


we  see  that  our  notion  of  solution  for  the  second  order  equations 
yields  p  and  B  L^{Q).  However,  the  conditions  inherent  in  (135)  when 
appropriately  decoupled  require  more  of  p  and  B.  For  example,  consider 
the  condition, 


-Vp+VxBxBq  =  f  e  L^W  (159) 


Take  the  curl  of  both  sides  in  D'(r2),  to  get 
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V  X  (VxBxBq)  =  V  X  f  e  H"k^)  (160) 

From  the  vector  identity  (valid  in  D'), 

V  X  (axb)  =  b«Va  -  a«Vb  -  bV«a  +  aV«b  ,  (161) 

we  deduce, 

B«V(VxB)e  H-1(J2)  .  (162) 

Also,  take  the  dot  product  of  (159)  with  B^,  to  obtain, 

Bq.Vp  =  Bg'f  SL^  (a)  (163) 

The  uniqueness  of  solutions  to  (89-90)  guarantees  that  there  will 

exist   gj^(x)  and  g2(x)  (which  amounts  to  choosing  an  appropriate  subset 

of  initial  data)   for  which   (159)   and   (163)   hold.   A  direct   not 

functionally  dependent  definition  of  the  appropriate  g^^'s  and  g2's  does 
not  appear  to  be  available,   however.    As   is   clear   from   the  above 

discussion,   however,  f  =  g  =  0  will   not  ensure  that  (159)  and  (163) 
hold. 
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